Length control of microtubules by depolymerizing motor proteins 
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In many intracellular processes, the length distribution of microtubules is controlled by depoly- 
merizing motor proteins. Experiments have shown that, following non-specific binding to the surface 
of a microtubule, depolymerizers are transported to the microtubule tip(s) by diffusion or directed 
walk and, then, depolymerize the microtubule from the tip(s) after accumulating there. We de- 
velop a quantitative model to study the depolymerizing action of such a generic motor protein, 
and its possible effects on the length distribution of microtubules. We show that, when the mo- 
tor protein concentration in solution exceeds a critical value, a steady state is reached where the 
length distribution is, in general, non-monotonic with a single peak. However, for highly processive 
motors and large motor densities, this distribution effectively becomes an exponential decay. Our 
findings suggest that such motor proteins may be selectively used by the cell to ensure precise con- 
trol of MT lengths. The model is also used to analyze experimental observations of motor-induced 
depolymerization. 

PACS numbers: 05.40.-a; 87.16.-b; 87.16.Nn 



I. INTRODUCTION 

In eukaryotic cells, microtubules (MT) are one type of 
cytoskcletal filaments, which perform several roles: they 
serve as tracks for intracellular molecular motor trans- 
port, provide structural rigidity to the cell and assemble 
to form a spindle during metaphase for the purpose of 
separation of duplicated chromosomes to the mother and 
daughter cells. MT are highly dynamic. In in vitro sit- 
uations, it is observed that a growing MT can suddenly 
start shrinking; this "catastrophe" is triggered by the 
loss of the GTP cap by hydrolysis of GTP molecules at- 
tached to the tubulin subunits of MT. A shrinking MT 
occasionally get "rescued" and start polymerizing again. 
This unusual process of polymerization and depolymer- 
ization is referred to as 'dynamic instability' [lj. 

The situation inside the living cell is much more com- 
plex and the MT kinetics in vivo is regulated by many 
proteins. In particular, several motor proteins belonging 
to the kinesin family are now known to function as MT 
depolymerizers 0, H, 0; IE] and are crucial in the forma- 
tion and maintenance of the mitotic spindle [f|. These 
include XKCM1/MCAK (mitotic centromere-associated 
kinesin) and Kif2A belonging to the kinesin-13 fam- 
ily and Kip3 proteins belonging to kinesin-8 family. 
MCAK is known to be active at kinetochores (the pro- 
tein structure which facilitates MT attachment to chro- 
mosomes) , whereas Kif2A is associated with centrosomes 
and Kip3p regulates microtubule-cortical interactions. 
Depletion/inhibition of XKCM1/MCAK has been shown 
to affect spindle length and the poleward motion of chro- 
mosomes during anaphase, while deletion of kincsin-8 
proteins leads to defects in positioning of the spindle 0. 
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However, the mechanism of depolymerization of MT by 
these motor proteins is only incompletely understood [7J . 

Recent in vitro depolymerization experiments with 
surface-immobilized MT have demonstrated that the de- 
polymerizing kinesins MCAK Q and Kip3p @, fiojuse 
a 'reduction of dimensionality' mechanism to target 
the MT tips for depolymerization: after binding non- 
spccifically to the MT surface, these proteins use dif- 
fusion (MCAK) or directed walk (Kip3p) to reach the 
tip(s), and their accumulation induces depolymerization 
of the MT from the tip. The depolymerization rate is, 
in general, length and motor-concentration dependent; in 
addition, MCAK and Kip3p are found to have vastly dif- 
ferent residence times on the MT. These results indicate 
that depolymerizing kinesins could be used by the cell 
for precise length regulation of microtubules. 

The central question of interest for us here is: what 
is the nature of the length-distribution of a set of micro- 
tubules in steady state, in a solution of free tubulin and 
depolymerizing motor proteins (henceforth referred to as 
depolymerizers for brevity)? In particular, it is impor- 
tant to understand how the depolymerizers contribute to 
the formation of a highly ordered structure such as the 
metaphase spindle, where most of the MT would need 
to be close to the spindle length. This aspect becomes 
significant, if we recall that the steady state length distri- 
bution of a set of microtubules undergoing dynamic insta- 
bility under standard in vitro conditions has an exponen- 
tially decaying form, without any preferred length [111 ]. 

In this letter, we present a theoretical framework for 
understanding the depolymerizing activity of a generic 
motor protein as described above. After deriving expres- 
sions for the concentration profile of bound motors on 
a MT, we calculate the rate of absorption of motors at 
the MT tips. We then construct the rate equations for 
the length distribution(s) of MT, which are then solved 
perturbatively in the limits of low and high motor con- 
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FIG. 1: An illustration of our model, showing the important 
kinetic processes. Depolymerizer motor proteins attach to 
the surface of a MT and are transported along the MT by 
diffusion and/or directed walk. The motors are trapped at 
the MT tips, where they induce depolymerization of the MT. 
The MT grows at a motor-free tip by adding tubulin subunits. 
Bound motors may detach from the surface or the tips of the 
MT. 



TABLE I: A list of experimental parameter values for 
MCAK [8|] and Kip3p [Sjj depolymerizers. The depolymer- 
ization rate is a function of length and motor concentration. 



Quantity 


MCAK 


Kip3p 


D 


0.38 /im's -1 




Vo 




3.6 fim min - 


Vd 


< 4/nm min -1 


< 2/im min - 


~k+ 


0.64nM - Vm -1 s -1 


assumed same 


k- 


1.21 s -1 


~ 0.004s -1 


K 


0.5s" 1 


0.03 s -1 



centrations. We also show that the model can be used 
to analyze available experimental data on motor-induced 
depolymerization. 



II. MODEL DESCRIPTION 

Since we are interested in length-related properties, we 
consider here MT of finite length (unlike earlier theoret- 
ical approaches [1, [II] to the problem, which concen- 
trated on semi-infinite MT). In our model, a depolymer- 
izer binds to a MT at a rate fc+co per unit length (co 
being the depolymerizer concentration in solution), and 
detaches at a rate fc_. A bound generic depolymerizer 
undergoes a biased random walk with drift speed Vq and 
diffusion coefficient D. Once a depolymerizer reaches a 
tip of the MT, it depolymerizes the filament at a rate v s 
(measured in /im min -1 ) while remaining attached to the 
tip. The depolymerizers accumulate at the tip, and de- 
tach from there at a rate k, which characterizes its proces- 
sivity(the tendency to stay attached to the MT tip while 
depolymerizing it). The tips, therefore, act as sinks for 
the depolymerizers which could accumulate there. The 
net rate of depolymerization, Vd, is a stochastic quantity, 
and depends on the number of accumulated motors. We 
further assume that a GTP cap is always present at a 
MT tip so that a depolymerizer-free tip will grow at a 
rate v g by adding tubulin sub-units (i.e., GTP hydrolysis 
is assumed to be sufficiently slow) . An illustration of the 
model is provided in Fig. [TJ 

The MT-bound depolymerizer density profile will be 
denoted by c(z, t) with < z < i , £ being the length of 
the MT at time t. The kinetics of c(z,t), in a frame of 
reference attached to the moving MT tip is given by the 
equation, 



dc „<9 2 c dc 



Ot 



dz 



(1) 



where v = v — v. m and v m — v g for a growing MT and 
v m = —Vd for a shrinking MT. Note that steric exclu- 
sion effects between motors have been neglected here for 



simplicity. In general, Vq depends on c; however, at suf- 
ficiently small c, modification of the last term of eq. ([TJ) 
to incorporate steric exclusion is equivalent to addition 
of just a constant to the detachment rate fc_(l2|. 

So far as the surface transport on the MT is con- 
cerned, the MT tips act as "traps" for the depoly- 
merases; however, a depolymerase can escape from the 
"trap" only by detachment from the MT during its MT- 
depolymerizing activity. Therefore, we treat the ki- 
netics of the tip-absorbed motor population separately 
from the MT surface-bound (mobile) density of motors. 
The boundary conditions on c(z,i), which are needed 
to solve eq. JT]), should reflect the "trapping" of motors 
at the tips, and the simplest case is absorbing bound- 
ary conditions: c(0, t) — c(£, t) = for all t. Us- 
ing time-independent boundary conditions requires that 
length fluctuations due to growth/shrinkage are suffi- 
ciently small, and the conditions for the same are derived 
now. We consider the v — case first. A single motor- 
typically spends a residence time r r ~ l/fc_ on a suffi- 
ciently long MT (I ^> \ADr r ), or otherwise gets absorbed 
at a tip within a shorter time interval ~ £ 2 /D. Clearly, 
only those motors that bind to the MT within a 'deple- 
tion zone' of length id ~ y/Dr r from a tip get absorbed at 
the tip, within a time interval ~ ^/-D ~ r r , the rest will 
get detached before they can reach the tip. The length 
change over r r is 5t ~ v m T r with v m = v g or — v s , and 
the condition \5£\ <C I is satisfied by \v m \ -C WDkZ. In 
the case of a motor which undergoes directed walk with 
velocity vo, the corresponding condition turns out to be 
\v m \ < vo- 

Table 1 lists the various experimentally measured pa- 
rameters for the two depolymerizers MCAK and Kip3p. 
For the purely diffusing MCAK, the above condition is 
satisfied for \v m \ <t^. 17\ra\ min -1 , and the analysis based 
on fixed boundary conditions should work well for growth 
rates of physiological interest. For walking Kip3p, the 
condition is marginally satisfied for low growth/shrinkage 
rates of \v m \ <C lpm min -1 . Nevertheless, we carry our 
analysis for the general case of a depolymerizer with di- 
rected motion as well as diffusion, assuming that the con- 
dition of small length fluctuations is satisfied. In partic- 
ular, we will henceforth assume that v ~ vq in eq. JTJ). 
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III. MOTOR DENSITY PROFILE AND 
ABSORPTION RATES 

The steady-state density profile c(z) of the motors, ob- 
tained from eq. ([T]), under the absorbing boundary con- 
ditions is 



c{z) = K 



1 



1 



smk{(3vl/2D) 



v(e-z) 

e 2D sinh(/3wz/2D) 



sinh(f3v(£ - z)/2D) 



(2) 



where = \/l + 4Dk-/v 2 and K = fc+co / fc_ . We em- 
phasize that the function c(z), being the density profile 
of only the motors not absorbed at the MT tips, vanish 
at z — and z = £■ However, the total concentration 
of the depolymerizers at the MT tips is non-vanishing 
because of the accumulation of the absorbed depolymer- 
izers there. We also note that eq. is invariant under 
the combined transformations v — > —v and z — > £ — z. 

It is interesting to look at the limiting behavior of the 
density profile in eq. @ in the limits v — > 0, D > 
(pure diffusion) and D — > 0,v > (pure walk), since 
these special cases are of experimental interest. In the 
first case, the distribution has the following well-defined 
form, which is symmetric about z = £/2: 



Cd(z) = lim c(z) = K 



1 



smh(£/X d ) 



sinh(z/A rf ) 



sinh((i - z)/X d ) 



; A d = y/D/k. (3) 



In the second case, after taking the limit D — > 0, we 
find that 



c v (z) 



lim c(z) 



K[l-exp(-z/X w )] ;X V 



v/k- 
(4) 

Note that eq. (gj) is non- vanishing expected 
for a pure directed walk towards the plus end. This does 
not contradict the imposed boundary conditions, because 
we still have lim^^o lim x >z c(z) = 0. 

Since the kinetics at the minus end of a MT is typically 
much slower than that at the plus end, we concentrate 
only on the plus-end kinetics for the rest of this paper, 
though motor accumulation is allowed to occur at both 
ends. The rate of absorption of motors at the plus-end 
is given by v(t) = —D(dc/dz) z= £, and has the general 
form: 



vK(l + f3)/2. 



(6) 



For pure diffusion and pure walk, eq. ^ reduces, re- 
spectively, to the limiting forms, 



and 



v d (£) = ~k+c X d tanh(£/2A d ) ; v -> (7) 

= k + c X lo [1 - exp(-£/X w )] = vc v (l) \D — > 0. 
MT {£ < A d or A,, 



(8) 

For short MT (£ < X d or X w respectively), v{£) ps 
^k + £ in the first case and v(£) w k + £ in the second; the 
difference of the factor of 2 arises from the fact that a 
diffusing depolymerizer can target either of the tips. In 
the opposite limit of long MT, both the rates approach 
their saturation values k+c$Xd and fc+coAm, respectively, 
(which may also be derived directly from eq. (J6j) ) . 



IV. RATE EQUATIONS 

We will now study how a given depolymerizer would 
affect the length distribution of a set of MT, when 
a steady state is reached by a balance between de- 
polymerizer binding/detachment and MT polymeriza- 
tion/depolymerization processes. For simplicity, we ne- 
glect the three-dimensional structure of the MT filament 
and imagine the MT as a linear polymer, made of sub- 
units of length b. We denote by P n (m,t) the fraction of 
polymers with m sub-units and n absorbed depolymeriz- 
ers at the plus-end. Let p g = v g /b be the probability per 
unit time for attachment of a subunit to a free tip and 
p s = v s /b be the rate of removal of subunits per motor 
(we assume henceforth that the rate of sub-unit removal 
increases linearly with the number of motors). The rate 
equations for P n {m,t) are as follows: For m = 1, 



dP (l) 
dt 

dP 1 (l] 



dt 



-p g P (l)+KPx(l), 

= P.Pi(2)-kP 1 (1) + i/(&)P (1), 



(9) 



while for m > 2 and n > 1, 



v{£) 



2wah{flvl/2D) 



P[ cosh(/3v£/2d) - e vl/ 



2D 



+ 



$mh{pv£/2D) 



.(5) 



v{£) vanishes at £ = 0, and is a monotonically increas- 
ing function of £ for all values of v and D. The saturation 
value for the general v > 0, D > case is given by 



dPg(m) 
dt 

dP n (m) 
dt 



= -p g [P (m)-P (m-l)} + 
nPi(m) — is(mb)P (m) 

= np s [P n (m + 1) - P„{m)} + 



(n + l)nP n+1 (m) + ^(TO&)F„_i(m) 
(tik + v{mb))P n (m) 



(10) 
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We now focus on the steady state of the model. Since 
a general analytic solution is difficult, we will adopt a 
perturbative approach, and look at the solutions in the 
limit of low and high depolymerizer concentrations. The 
transition between these regimes is controlled by the di- 
mensionless ratio 77 = i/ max /K which may be used as the 
small parameter in the low density perturbation expan- 
sion. This is seen through the following argument. 

Let us define a distribution of the number of accu- 
mulated motors: \ n = J2m=i Pn( m ) ■ It is convenient 
to use the numbers a n < 1, defined through the re- 
lation Emd^W^f™) = fmaxCKnXn- From eq.( 02) 
and eq.( [TU)) . it can be shown that, in steady state, 
Xn follows a Poisson-like distribution of the form \n — 
Xo(v n / n *) n"=o ai f° r n — 1' wnere Xo is determined 
through normalization: Y^n°=o = L 

Low motor concentration: For r\ <C 1, we can neglect 
P„ with n > 2 in eq. (|10|) , in comparison with Pq and 
Pi. By keeping terms upto 0(rj), in the continuum limit 
(to 3> 1, £ ~ mb), we arrive at a single combined equation 
for Pq and P\ in steady state: 
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FIG. 2: The distribution in eq. (I17|l (normalized such that 
J Po(£)d£ — 1) is plotted against the dimensionless 'length' 
x = ^/2A,j(eq. Q) for three different values of ( = r\jr\ c . 
The other parameters were fixed at the MCAK values (table 
1). We assumed v s — Vd and v g ~ lfim min -1 . From table 1, 
Xd ~ 0.56/im for MCAK. Inset: the mean length as a function 
of ( (computed numerically, 8x = 10~ everywhere). Using 
the parameters in table 1, the critical MCAK concentration 
corresponding to 77c is 1.39n 9 /« s nM. 



dP 



dP x 



which has the general solution (for £ ^> b) 



(11) 



PM) = Cexp (—1-— [ v(£')dA = —Pi{£) (12) 
\v s v g J J Vg 

where C is a constant of normalization. In the absence 
of a boundary, this solution is normalizable only when 
77 > ?7 C , where rj c = v g /v s implicitly gives a critical motor 
concentration, below which a well-defined steady state 
distribution does not exist. The above result for r; c is, 
however, strictly true only when v g /v s <C 1 since the 
present analysis assumes rj <C 1. In the more general 
case, a steady state should still exist at sufficiently large 
to but the critical concentration would be given by a 



f{v g /v s ), with f(x) 



x as 



relation of the form r\ c 
x 0. 

High motor concentration: To analyze the case rj 3> 1, 
it is convenient to re-express eq. © and eq. (fTTJj) (in the 
continuum limit) as 



dP 



(0) 



= f(£){Pi°l 1 (£)-Pi°\£)] (14) 



Of 



and P (0) W = P^iOjexpi-Vg 1 f£d£'v{£')) from 
eq. (fT5)) . The solution of eq. (|14[) is found after a few 
elementary calculations: 



where 



p£\l)=p£)(Q)e-^ti M ' v{t '\ 



p(°)(0) = P W(0)n 



i Vg + mv s 



(15) 



(16) 



The unknown constant P^ is determined by normal- 
ization. Eq. (|15[) is a monotonically decreasing function 
of £, however, non-monotonic terms appear in the first 
order and above. 



v g 9P 

An ax 3^ 



r)- l Pi{Z)-f(l)m) 



nv s dP n 

^max &£ 



(n + l)n- l P n+1 {i) + /(Z)P„_i(*) - 

(nr,' 1 + f(£))P n (£),(13) 



where the function f(£) = v{£)/v ma ^ < 1. We now 
expand all P n (l) in a perturbation series in ?y _1 : P n (£) = 

J2m=o V~ m Pn m ^ {£) for n > 0, and then substitute into 
eq. (I13[) . The zeroth order distributions are then found 
to satisfy the iterative equation, 



V. APPLICATION TO SPECIFIC CASES 

Our results so far have been general, within the lim- 
its of validity of the assumptions stated in the second 
section. We will now examine their implications for the 
specific depolymerizer motor proteins studied in experi- 
ments. 

The solution obtained in cq. (Ti"2"|) in the limit of low 
concentrations is a non-monotonic function of the length, 
increasing exponentially as ~ e ( K -/ v ^) i for small £, and de- 
creasing exponentially at large £, with a peak at an inter- 
mediate value. The location of this peak depends on the 
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FIG. 3: Similar to the previous figure, using the distribution 
in eq. psp and Kip3p parameter values (table 1). Here, x = 
l/\ w (eq.(|4]))with \ w ~ 15/^m for Kip3p. The critical Kip3p 
concentration corresponding to r\ c is ~ 0.003v g /v s nM. 
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FIG. 4: The ratio of the standard deviation to the mean 
length of the distributions in figO and fig0] as a function of 

C = v/vc- 



ratio k/v s , and for fixed v s , the peak is more pronounced 
for large k (low processivity) and vice- versa. The depoly- 
merizers MCAK and Kip3p have very different processiv- 
ities (table I), and it is therefore interesting to look at 
these specific cases in more detail. Eq. (fT2|) reduces to 
the following forms in the pure diffusion (v — > 0, D > 0) 
and pure walk (v > 0, D — > 0) limits respectively: 



P (£) oc e Rl ' v ' [co8h(^/2Ad)] 



Po{i) oc exp 



(«o = 0) (17) 



-t/x w 



(D = 0) 
(18) 

Fig. [2] and fig. [3] show the normalized forms of the dis- 
tributions in eq. (fT7|) and eq. (fl8)) respectively, plotted 
against the dimensionless length variables x = £/Xd and 
x — £/X w , for three values of C, = rj/rj c > 1 (the empirical 
values of k,D,vq and v s are taken from Table 1). Al- 
though, as we explained in the beginning, our theory is 
more applicable to MCAK than Kip3p, the difference be- 
tween the two cases is nevertheless striking: the MCAK 
distribution shows a peak, while the Kip3p distribution 
is essentially a monotonically decreasing function (the 
peak is too close to the origin to be visible in the plot). 
The contrasting behavior is primarily due to their very 
different processivities. 
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FIG. 5: The figure shows the mean number of tip-absorbed 
motors (computed numerically using eq. I19|l . plotted as a 
function of the MT length for MCAK, for two values of the 
bulk motor concentration cq. The inset shows the saturation 
value (computed at t — 5/im) plotted against Co. 



In order to characterize the MCAK and Kip3p-induced 
distributions better, we also looked at how the lengths are 
spread about the mean value in each of these cases. In 
fig.[U we plot the relative MT length fluctuation Ax/ (x) 
for MCAK and Kip3p, as a fu nction of th e dimensionless 
ratio C = v/Vci where Ax — W (x 2 ) — (a;) 2 is the standard 
deviation. The relative fluctuation for MCAK is typically 
smaller than I, and, interestingly, also shows a minimum 
of ~ 0.3 at £ ~ 1.37. For Kip3p, on the other hand, the 
relative fluctuation increases rapidly with £ and saturates 
at unity, which is characteristic of a purely exponential 
distribution. MCAK appears to produce a tighter control 
of length than Kip3p because of its lower processivity. 



VI. LENGTH-DEPENDENT 
DEPOLYMERIZATION 

We will now focus on the analysis of the existing ex- 
perimental results on motor-induced depolymerization 
within our model. The rate of depolymerization of a 
microtubule depends on the number n of motor proteins 
accumulated at the tip. Experiments with fixed micro- 
tubule length measure the mean depolymerization rate 
(Vd) = v s {ri)£, where (n)e gives the mean number of mo- 
tors attached to a microtubule of length I. To calcu- 
late this average, we define the distribution Qi(n) for the 
number of motors attached to the tip of a MT of length 
I. We assume that a motor protein initially bound to 
one protofilament will continue to travel along it until it 
reaches a tip, without hopping between filaments. In this 
case, absorption at a certain protofilament tip is possible 
only if that particular tip is free, the probability of which 
is given by (I — n/n*) 1 where n is the total number of 
motors attached at the MT tip and n* is the maximum 
number of motors that can be absorbed at a given time. 
Since a MT usually has 13 protofilaments, we assume 
n* = 13 for concreteness. 

The rate equations for Qe(n) are as follows. For motor- 
free tips, we have 
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FIG. 6: The average number of tip-absorbed motors in the 
presence of Kip3p, for three different values of Co- Note the 
difference in length scales between this and the previous fig- 



We now turn to the case of Kip3p, in which case a 
strong length-dependence of depolymerization was ob- 
served in experiments @ . Fig. [B] shows the theoretical 
plot of {n)i against I for three motor concentrations. The 
length-dependence here is much stronger than MCAK, 
and for low concentrations, saturation is not reached 
even for i = 20/im. These observations qualitatively 
agree with the corresponding experimental results. How- 
ever, a close quantitative agreement (in particular, ex- 
periments show a sharp rise in depolymerization rate be- 
tween co=3.3nM and co=5.8nM) is missing in this case. 



VII. CONCLUSIONS 



dQe(0) 
dt 

dQejn) _ 
dt 



1)W 



= KQ t {l) - u(£)Qt(0) 
■■ (n + l)KQi(n + 1) + 
Qi(n — 1) — Qi(n) + 



Qi(n — 1) — nnQi(n) ; 1 < n < n* 

dQejn*) m is(£) * 

— = -n nQi(n ) H —Qi\n - 1) (19) 

at n* 

The steady state solution of eq. (flU)) is given by 



Qt(n) = Q e (0)- 



1 (v{l) 



-) l<n<n* 



(20) 

where Qe{0) is fixed by normalization: J2n=i Qt{ n ) = 
1. The mean number of tip-accumulated motors is then 
given by (n) e = YZ=i nQe(n). 

In fig. O we have plotted numerically computed (n)( 
as a function of £ in the case of MCAK for two differ- 
ent motor concentrations cq. It is found that apprecia- 
ble length-dependence of {n)t, and hence the mean de- 
polymerization velocity (vd), is found only for t < 2/zm. 
This is in agreement with experiments, where no length- 
dependence of depolymerization was observed in the case 
of MCAK, for MT longer than 2//m. The saturation 
value of the depolymerization rate (plotted in the inset, 
for I = 5/Ltm) depends strongly on the motor concentra- 
tion, and this curve also agrees well with experimental 
results Q and previous theoretical predictions [r2| . 



To conclude, in this letter, we have formulated a gen- 
eral theory for the action of MT-depolymerizing motor 
proteins. Specifically, we considered the limit where the 
motor density profile becomes stationary much before 
the MT length distribution reaches its steady-state. We 
showed that the rate of accumulation of the motors at the 
MT tips, and consequently their depolymerizing activity 
itself, is strongly length-dependent. This has a rather 
pronounced effect on the length distribution of the MT, 
which displays a peak before decaying exponentially at 
large I. Interestingly, the processivity of the depolymer- 
izer plays an important role in determining the nature 
of the length distribution: a depolymerizer with low pro- 
cessivity produces a more pronounced peak in the length 
distribution, and is likely to be more useful for precise 
length regulation. 

Our theory is relevant to future experimental studies 
on the role of depolymerizing motor proteins in length 
regulation of MT, especially in the context of formation of 
the metaphase spindle. Indeed, very recent experiments 
using fluroscent speckle microscopy have shown that the 
length distribution of individual MT in a meiotic spindle 
is strongly non-monotonic [l3j . It would be interesting if 
the predictions made in this paper could be put to test 
in future in vitro experiments with depolymerizing motor 
proteins, where the MT lengths are carefully monitored. 
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